C-axis Transport in Bilayer Cuprates and Relation to Pseudogap 
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We consider the effect of interband transitions on the c-axis conductivity, DC resistivity, and ther- 
mal conductivity in a plane-chain bilayer model of a cuprate. The relation between the c-axis 
resistivity and thermal conductivity is governed by the Wiedemann-Franz law. When the perpen- 
dicular hopping matrix element between chain and plane {t±) is small, the c-axis DC resistivity 
shows a characteristic upturn as the temperature is lowered and the infrared conductivity develops 
a pseudogap. As t± is increased, intraband transitions start to dominate and a more conventional 
response is obtained. Analytical results for a simple plane-plane bilayer are also given, including the 
frequency sum rule of the optical conductivity. 



PACS numbers: 74.80.Dm, 78.30.-j, 74.25. Jb, 74.72.-h 



I. INTRODUCTION 

In a previous paper ||l[, we derived expressions for, 
and computed, the normal-state optical conductivity for 
a system of two layers per unit cell coupled through a 
hopping matrix element t± which introduces the third 
dimension. In addition to the usual intraband terms, 
which give a Drude-like response, interband transitions 
are now possible and these have quite a different temper- 
ature and frequency dependence. It was found that for 
sufffciently small values of t± , the interband terms, which 
are of order i^ , dominate over the intraband contribution 
which goes like t^ for the c-axis. In this limit, the c-axis 
response can be very different from Drude. In general, 
the interband transitions occur at energies which extend 
over the entire band width and structure in the c-axis 
conductivity can be directly related to certain features 
of the band structure. For the particular case of a Cu02 
tetragonal plane and a set of CuO chains ^-^, it was 
found that the normal-state c-axis response can develop 
a gap or pseudogap at low lu depending on the details of 
the two Fermi surfaces involved. 

Estimated values of t± have been provided by Zha, 
Cooper and Pines ||]. They give tj_ ^ S.OmeV for under- 
doped YBCOg.es and 30 to 40meV for optimally doped 
YBCOe.Qs- These estimates are much larger than for 
the more anisotropic compounds such as Bi2Sr2CaCu208 
where t± is perhaps as small as O.lmeV. It is clear 
from these estimates that underdoped YBCO could be in 
the interband dominated regime and show a pseudogap 
response while optimally doped and overdoped YBCO 
would show a more usual Drude-like c-axis response as 
in this case, intraband transitions would dominate over 
interband. 

In this paper, we wish to consider the c-axis transport 
and study the effect of interband transitions on its tem- 
perature dependence. For fully oxygenated YBCO, the 
resistivity p{T) in the c-axis p-[l^ is found to exhibit 
metallic behavior and show a nearly linear temperature 
dependence which follows that in the a-b plane but, of 
course, has a very different magnitude with Pc/ Pab of or- 



der 30 to 70. On the other hand, for deoxygenated under- 
doped samples corresponding to small t±, the situation 
is very different and Pc(T) shows a semiconductivity-like 
behavior |p|,p^-p^ at low temperatures {dpc/dT < 0), 
while at high temperatures, the linear behavior is re- 
tained. The ratio pd Pab is however very much larger 
than that observed in the fully oxygenated case with 
Pc/ Pab &t T — lOOX of order 10'^. Estimates of c-axis 
mean free paths based on such large values of pc lead to 
values that are of the order and even less than the c-axis 
spacing between planes. This realization has led some to 
suggest that conventional three dimensional Bloch trans- 
port breaks down and that c-axis transport proceeds 
through an incoherent mechanism due to confinement 
coming from the highly two dimensional nature of the 
electronic states. Here we consider that the interband 
transitions begin to dominate when t±^ is small. This 
effect has its origin in the fact that YBCO has several 
atoms per unit cell. In the real situation, there are two 
copper-oxygen planes and one set of chains. For simplic- 
ity, here we will consider only a plane-chain model which 
is believed to capture some of the qualitative features of 
the interband transitions. 

In Sec. II, we establish the Kubo formula for the AC 
infrared conductivity in the case of a plane-chain model 
coupled through a transverse hopping matrix element tj_ 
in tight binding. The usual formula becomes a trace 
over the two bands and the electromagnetic vertices are 
now a 2 X 2 matrix with off-diagonal elements related 
to the derivatives of the Hamiltonian matrix. The re- 
sulting formula easily separates into interband and intra- 
band parts. In Sec. Ill, analytic results for the infrared 
conductivity and the DC resistivity are presented for the 
case of infinite bands with minimum energy difference 
between them of A. One band is taken to have a one- 
dimensional character and the other two-dimensional. 
The role of the pseudogap A in these quantities is made 
clear. In Sec. IV, the thermal conductivity is consid- 
ered and the Wiedemann-Franz law is established for the 
c-axis including interband as well as intraband contribu- 
tions. In Sec. V, we treat the case of a plane-plane bilayer 



which is important because simpler analytic results can 
be obtained including for the conductivity sum rule. In 
Sec. VI, tight-binding bands are introduced and numeri- 
cal results presented for this more complicated case which 
does not yield to any analytic analysis but includes such 
effects as the Van Hove singularity in the electronic den- 
sity of states. Nevertheless the results show no qualita- 
tive difference with those presented in previous sections 
based on simplified infinite bands. A short conclusion 
can be found in Sec. VII. 



II. FORMALISM 

We start with the Hamiltonian for a plane-chain bi- 
layer coupled through a transverse matrix element t±. 
In a matrix notation, we start with creation C^ and an- 
nihilation Ck operators 






cl- 



a = 



Clk 
C2k 



(1) 



with c{[(;(''2k) creating an electron in the state of mo- 
mentum k in layer 1 (2), respectively. The Hamiltonian 
operator Hg is 



Ho^Y. ^k^(k)Ck, 



where the Hamiltonian matrix 



hik) 



t{K) 6(k) 



with 



a(k) 

6(k) 






2m 



k^- fJ; 



for the plane and chain respectively and 



t[k^) = -2t_LC0S ( — 



(2) 

(3) 

(4a) 
(4b) 

(5) 



with d the distance between chain and plane in the z- 
direction and t± the couphng strength. In (^), fj. is the 
chemical potential, A is the minimum energy difference 
between the two bands and is responsible for the pseu- 
dogap in the c-axis response, H is the Planck's constant 
over 27r and m is the bare electron mass. 

When a vector potential A(c|) is applied to the system, 
the Hamiltonian changes to Hi^q] 



H^Hn- 



2mc 



pA 



(6) 



to the lowest order in A, with the vector 



fE^i 



dh(k) 
dk 



a 



(7) 



Using the Kubo formula for the current-current correla- 
tion function, this leads directly to a real-part conduc- 
tivity in the c-axis of the form (Reazzi'^) = Cc(^)) 



adio) 



o 2^ / 'oZ 



dx f{x) — f{x + huj) 



^ ^ .l-r^ 27r uj 



xTr[A(k, x)7,(k, k)A(k, X + nw)-i,{k, k)], (8) 

where "Tr" is a trace, A{k, uj) is the spectral function 
matrix, and 7z is the associated c-axis electromagnetic 
vertex function matrix given by 



1 dh{k) 



h dkz 



(9) 



The "trace" operator in (g) allows one to work in any 
frame which is convenient. We will choose a frame in 
which the Green's function (and hence the spectral func- 
tion) matrix is diagonal. Consequently, one obtains (fol- 
lowing Rcf. |§) 

n^ J ^2^ '^ 



crc(t^) = e^ 



(10) 



(AiiA'h + ^22^^2)7?! + (^11^22 + ^22^'ll)7?2 



where An = Aii(k, x) and A'^^ = ^^(k, x + hoj) and the 
vertices 

2t 

711 = vi_ 



712 = v± 






(11) 



with v± = 1/hdt/dkz — t±d/hsm(kzd/2). More explic- 
itly, in the normal state, the spectral functions 

2r 

A„(k,.).^— -^, (12) 

where ci = e+, e2 = c-, and Ti is the total scattering 
rate in band i. A realistic model for the high- Tc oxides 
is that the scattering rate is dominated by the inelastic 
scattering with T linear in temperature T. The energies 



e.^^-^±J(K^ 



t^ 



(13) 



are the two renormalized bands. 

It is clear from the form of Eq. (noh that with the initial 
and final states in the same band, the first term corre- 
sponds to the intraband contribution and the second term 
corresponds to the interband contribution with the ini- 
tial and final states in different bands. From the form of 
the vertices given in (^l|), one sees immediately that the 
interband contribution (^ t^) dominates the intraband 
contribution (~ t\) in the small-t^ limit. For simplicity, 
we will assume the scattering rate Fi = F2 = F in all our 
calculations. 



III. CONDUCTIVITY AND DC RESISTIVITY 



Substituting (y), (|ll|), and (|12|) into (|10|), the intra- 
band contribution for the c-axis conductivity is exphcitly 
given by 



^intra 



(u) 



^E«5 



n 



2t 



fix) - fix + htu) 



2r 



2r 



dx 

2^ 

2r 



(x-e+)2+r2 ix + huj-e+)^+T^ 

(M) 



(a;-e_)2+r2 (x + ?iw - e_)2 + r2 



To reduce (041) further, we first replace the momentum 
k-sum by an integral 



n 



E 



X "27 



A^(0)de, 



(15) 



where e can be either e+ or e_ and a constant "2D" den- 
sity of states per spin at Fermi level A^(0) = m/27r?i 
is assumed (in the small-t^ limit). Secondly, we replace 
the factor e+ — e_ in (|lj) by the (constant) band gap A. 
This approximation is quite good and has only a small 
effect on the final results. As a result of the above sim- 
plification, the kz dependence in the integrand comes in 
only through tikz) and vj_ikz) which can be easily car- 
ried out to get J^^ ^Vi_t^ = t\d/2h?. Consequently, 
(|l4|) reduces to 



\u) 



_ 8eHlNiO)d /t±_ 



de 



dx (16) 



/(x)-/(x + M 



LU (a;-e)2+r2 ix + Hlj - e^ +T^' 

where the two bands give the same contributions for the 
intraband part of the c-axis conductivity. Furthermore, 
with the aid of 



de 



(a;-e)2-Hr2 ix + htu - e)^ + T^ 

2r 



and 



dx 



ihujY + 4r2 



fix)-fix + nuj) 



one immediately obtains the final result 
d>eH\NiO)d (t^ " 



\^) = 



2r 



h \AJ ihLjf + 4r2 

which is nothing but the "Drude" theory. 



(17) 



(18) 



(19) 



In analogy to (|lj), the interband contribution of the 
c-axis conductivity is given by 

dx fix) — fix + Tiuj) 
, 27r 



'{^) 



1 



k 



UJ 



2r 



2r 



(a;-e+)2+r2 ix + Ulj - e^)"^ +V^ 

2V 2r 



+ 



(x-e_)2+r2 (x + ^tj - e+)2 + r2 



(20) 



where we note that the second term dominates when 
a; > (because e+ > e_). To calculate (|o|), the simplifi- 
cation (pSl) used for the calculations of intraband part is 
no longer valid because now the two bands come in to- 
gether and the correlation and difference between these 
two bands are certainly important. One therefore needs 
instead to replace the momentum k-sum in ( |20| ) by the 
integral 



-y 

k 



" dk^ 

^ 27r 



(ikx 

'2^ 



Q/Ky 

"27 



(21) 



To simplify the calculation, we will approximate e+ = ^i 
and e_ = ^2 in the small-tj^ limit (see ([l3[)). This simpli- 
fication was not done in Ref. M where the k-sum is car- 
ried out directly by numerical integration over the first 
Brillouin zone. Note that this means that they are ef- 
fectively including the Van Hove singularity in the elec- 
tronic density of states which we have taken as constant. 
The integration over kz can be carried out first to be 
I-i% T^^l = t\d/2n^. Moreover, with 



dk 



1 



"^ (A-?i^fc2/2TO)2+r2 
1 



(22) 



dk 



™^°^ in (|2^) can be reduced to (changing the variable 

y ~* 



de- 



eHlNiO)d 



u2 



\/27r 



X 



dx 



fix) - fix + hu) 



UJ 



de^ 



-^ Ve- +Ai 



(a; + ncc;-e_)2+r2 y( 



A)2+r2 



^ix - e_ - A)2 



(23) 



+ -r 



- r2 - (x - e_ - A) 
F 



(a;-e_)2+F2 ^(x -e^ + Tioj- A)2 + F2 



\J^/ix -e^+hcu- A)2Tr2 - (a; - e_ + huj - A) 



with N{0) = m/2Trh . A crucial point seen in ( p^ ) is 
that the chemical potential (or filling) plays a role in the 
interband transition, as expected. 
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FIG. 1. c-axis conductivity as a function of photon en- 
ergy for a plane-chain bilayer at different temperatures with 
a small t± = 2meV. The minimum energy difference between 
the two bands for a given k in the first 2D Brillouin zone is 
A = 20meV. 

Based on the intraband contribution given by the an- 
alytical result ( |l9[) and the interband contribution given 
by the numerical calculations of pq), we have plotted 
the c-axis conductivity in Fig. yj for various temperatures 
at fixed tj_ = 2meV and in Fig. for fixed temperature 
T = 30K at various t±. We have used the same pa- 
rameters: A — 20meV, /i = SOOmeV, and F is taken to 
be linear in temperature with F = 20meV at T = lOOA' 
throughout this paper unless specified. This is motivated 
by the fact that the inelastic scattering rate is observed 
to be large and of the order of Tc |g] at the superconduct- 
ing critical temperature (Tc) and varying almost linearly 
with temperature T. Also the evidence is that good sam- 
ples show little or no residual resistivity and are therefore 
in the pure limit (there is no residual elastic scattering) . 
One sees in Fig. |] that a pseudogap can develop at low 
temperatures for small t± (in which case the interband 
transitions dominate). At higher temperatures, the pseu- 
dogap is invisible due to the stronger inelastic damping 
effect. 

While in this paper we have made simplifications and 
approximations to the band structure, for example, we 



used the infinite-band approximation (15) with constant 
density of states evaluated at the Fermi energy and taken 
it out of the integral, this was not done in Ref. [Q where 
the full sum over the Brillouin zone which includes a Van 
Hove singularity in the electronic density of states was 
performed numerically and the effect of t± on the bands 
(see Eq. (Oh) was fully included. We observe however 
that there are no qualitative differences between the re- 
suhs of Fig. 1^ and those of Fig. 2 in Ref. 0. The resuhs 
of Ref. [y do show some structure at higher energy which 



reflects special features of the plane-chain tight binding 
bands used in that study but not in this work, but these 
are small and we can safely proceed with a discussion 
of the DC resistivity using the simplified infinite band 
structure of Eq. (0) which allows the required numerical 
work to be greatly simplified. 
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FIG. 2. c-axis conductivity for a plane-chain bilayer for 
different values of coupling t^ at T = 30K with minimum 
energy between the two bands equal to A = 20meV. 

Fig. shows that by increasing strength of ij_, the 
pseudogap can disappear. For clarity of presentation, all 
intensities in Fig. have been divided by i^ . This leads 
to a nearly constant interband contribution plus a intra- 
band (Drude-like) contribution proportional to t'^ which 
increases with increasing values of t±. In the case of 
larger t± , the intraband transition and hence the Drude- 
like behavior becomes dominant (long dashed curve with 
t± = 20meV) over the smaller interband contribution 
which now appears above ~ 20meV. We stress and it is 
clear from the figures that the value of the real part of 
the conductivity at higher energies is entirely due to in- 
terband contributions since in our model the Drude con- 
tribution is sharply peaked about cj = on the energy 
scale of 150 to 200 meV. 

We have also plotted in Figs. ^ and || the c-axis DC 
resistivity for the plane-chain bilayer normalized to its 
value at 300i^. As shown in Fig. pi when tj_ is small 
the c-axis DC resistivity shows a characteristic upturn 
as the temperature is lowered such as those observed in 
semiconductors. This characterizes the importance of the 
interband transition. A "pseudogap" in the joint density 
of states (interband contribution) suppresses the inter- 
layer hopping and gives rise to the semiconductor-like 
Pc{T). When t±^ is increased, intraband transitions start 
to dominate and a more conventional response (linear in 
T) is obtained. A similar plot as given in Fig. ^ given 
in Fig. H shows that when the pseudogap is reduced (or 
relatively the t±^ is larger), the upturn feature becomes 
much less dramatic. 
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FIG. 3. Temperature-dependent c-axis DC resistivity for a 
plane-chain bilayer for different values of t±. The minimum 
energy between the two bands is A = 20meV. 

To gain more insight into the relationship between the 
c-axis conductivity and DC resistivity in a plane-chain 
bilayer, we derive in the following some analytic results 
under certain limits. When temperature is not too high 
and hence F is sufficiently small (it is assumed to vary 
linearly with T), one can approximate 



(x-e_)2+r2 



7r(5(e_ — x) 



(24) 



in the second term of (23) which one recalls dominates 
in the case of a; > which is of interest. Moreover, 
due to the behavior of the Fermi distribution functions, 
the integration over x in (|2^) contributes mainly from 
X — —huj — ksT to a; = ksT. One thus can replace 
/_ de- by J_ de^ assuming huj ^ /i. Consequently, 
after the integration over de- 
reduces to 



the second term of 122 



eHlN{0)d 



V2tt' 



dx 



f{x)-f{x + hiu) 1 



Vx+Jl 



sjihuj - A)2 + T^J^[hu: - A)2 + F2 - {hoj - A) 



(25) 



Eq. ( p5| ) allows one to discuss some analytic results of 
^mtcr ^g^j^j hence CTc)- 

For high frequencies [fiuj — A ::^ F, ksT), one can ap- 
proximate /_^ dx [f{x) — f{x + huj)] by J_f^^ dx, and, as 
a result, the associated x-integration can be performed 
easily to be 



dx 



1 



Next by expanding 



= 2(^- V/i-w). 



y/{hLU ~ A)2 -I- F2 - {hiJ - A) 



y/2{tlUJ - A) 



(26) 



, (27) 



one obtains following fe 



winter 



r 



^- ^[i~UJ 



(28) 



u}\/'hijj — A 

which is independent of F (or T) and is proportional to 
a;~2 when /z ^ w ^ A. 
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FIG. 4. Similar plot as Fig. |[ but for A = lOmeV. 

For low frequencies {Tllo <^ A), one has [f{x) — f{x + 
hLu)]/Lu ~* h6{x), and following ([2^) 



'H 



tlfi-^- 



(29) 



^/{hio - A)2 + T^^./{hL7- A)2 + F2 -~{n^- 



A) 



which is proportional to F when F <C A and is propor- 
tional to F^2 when F ^ A. From the result (|29| ) and 
( p9| ) at (jj = 0, one can also get an expression for the DC 
resistivity {pc{T)) for the plane-chain bilayer, which is 
then given by 



PciT) 



ti^r 



(30) 



V8^VA2+F2v/\/A2Tr2 + A 

It is clear from (|3^) that for t± <$:. A the first term in 
the square bracket is small and becomes important only 
when the second one is also small. For F small, 



PciT) 



iT^r 



F2 



vVA^ 



(31) 



the Drude contribution to the c-axis resistivity Pc{T) goes 
like F while the interband contribution on its own has the 
inverse dependence and goes like F"-'^. This is the physics 
underlying the upturn in the resistivity of Figs. and 
as T is lowered and is also the reason the pseudogap 
forms as seen in Fig. n^. The two effects are intimately 
connected in our theory. 
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FIG. 5. Temperature-dependent c-axis thermal conductiv- 
ity for a plane-chain bilayer for different values of t±. 



KciT) ^ tl 



t± 



(34) 



r-^ 



^ ^ ' %/8 ^VA2 + rV\/A2 + r2 + A, 

and for small T (using (pi])) 



K,(r) ^ t2 



4 y^ 



(35) 



We have plotted Kc{T) normalized to its value at 30QK 
using (p4) in Fig. 0. For small tj_, the significant de- 
pression of Kc{T) at low temperature characterizes the 
importance of the interband transition in which case one 
will also see the depression of optical conductivity and 
the upturn of the DC resistivity associated with pseudo- 
gaps. 



V. PLANE-PLANE BILAYER 



IV. THERMAL CONDUCTIVITY 

In analogy to (ph for the c-axis optical conductivity, 
one can derive the c-axis thermal conductivity in terms 
of the Kubo formalism p6|,p7[| 






i d fix) 
dx 



cTr[A(k,x)7.(k,k)A(k,x)7,(k,k)], 



(32) 



calculated at a; — 0. Here f{x) is the Fermi distribu- 
tion function, A(k, x) is the spectral function matrix, and 
72 (k, k) is exactly the vertex function matrix for calcu- 
lating the current-current correlation function given in 
(0) . As we have done in ( pO| ) for the optical conductivity, 
due to the trace operator, one can work in a frame in 
which the Green's function is diagonal. This will lead to 
the fact that the c-axis thermal conductivity can be sep- 
arated into intraband and interband contributions. How- 
ever, this is not needed since we can establish that the 
c-axis thermal conductivity obeys the Wiedemann-Franz 
law. Due to the fact that the main contribution to inte- 
gration over X for both (^) and (||) at w = comes from 
the small region around x = Q, the 2:-dependence inside 
the trace for both (^2|) and (||) is removed. One easily 
finds by comparison that the ratio of the c-axis thermal 
conductivity to the c-axis DC electrical conductivity is 



Kr 



crc{uo = 0) 

2 

T 



dx X 



,df{x) 



dx 



°° dx ^^(^) 



e 



dy 2y sech y ■ 



ks 
e 



T. 



dx 

(33) 



The above result is simply the Wiedemann-Franz law. 
Using dot) with (H), one obtains (recall F - T) 



Even simpler analytic results, which shed further light 
on the physics of our model can be obtained for a plane- 
plane bilayer in which the two uncoupled bands are given 

by (|) with ^2 now replaced by 6 = ^(^^ + ^y) - M- 
In this case, one can also have analytical results for the 
interband part of the c-axis conductivity 



'■(^) 



N{0)d 



2r 



{TiLo -f A)2 + 4F2 



2r 



{TiLj - A)2 -f 4F2 



(36) 



The result (BQ) for crj,"*^'' for the plane-plane bilayer leads 
to a peak feature at w = A but is otherwise similar to 
the results given by (Eq) (or more generally (23)) for a 
plane-chain bilayer. The simple analytic expression ( |3q ) 
thus helps gain physical insight into the c-axis transport 
properties. 



A. Conductivity 

The total c-axis conductivity given by the sum of (M) 
and M) will be 



a,(^) = <*-(a;)-Hai"'-(c.) 

2 



eH\N{Q)d 



X 



4r 



2tj 

^ ; [huoY + 4r2 
2r 2r 



(37) 



,(?iw-HA)2-h4r2 (/iw - A)2 -h 4F2 
Some features are: 

• a™^^'^ ex t\, crj,"*""^ ex t\, and the relative intensity 

.riiitra /_iiitcr 



C^" /C^c « tj^ 



• When w -^ and T (or T) is small, crj.'^*'''^ - T'^ ~ 
T-^ (Drude like), while a^"^'"' ~ r - T (non-Dmde 
like). 

• The high-frequency tail of adu;) goes as w^^ and 
decays very fast. 

• The peak at a; = A corresponds to the pseudogap 
and by changing the ratio of tj^ to A, the pseudogap 
can disappear. 



B. Frequency Sum Rule 



Because of 








r°° 4r f°° 


2r 


2r 

1 


TT 


' (nw-A)2 + 4r2 


n' 


one can easily obtain 


r, , , eHlN{0)dn 

/ duj (Tciuj) = i-2 

Jo n 


-(^)1 


1 



(38) 



(39) 



which states the w-sum rule for the optical conductivity. 



C. DC Resistivity 

Using the result ( |37| ) in the limit of a; ^ 0, we obtain 
the c-axis DC resistivity 



VI. CASE OF TIGHT-BINDING BANDS 

So far we have used the simplified band structures (^) 
and have tried to get simple analytic final results. These 
have required some approximations which we now wish to 
remove. We start with tight-binding bands with momen- 
tum k confined to the first Brillouin zone of the copper- 
oxide plane. We take 



and 



^1 = —2ti[cos{kj;a) + cos{kya) 
—2B cos{kxa) cos{kya)] — ^i 



-2t2 cos(fcj,a) — fi2- 



(42a) 



(42b) 



For a model with {ti,t2,fJ.i,fJ.2,t±} ~ {70,100,-65, 
— 175, 20} in meV and B — 0.45, the corresponding Fermi 
surface contours are shown in Fig. O. The dashed line are 
the unperturbed Fermi surface which have been chosen 
so that they do not cross. They are also the contours 
for the hybridized case with finite t± and fcz = ^ in 
which case t{kz) in Eq. (ph is zero. The solid lines are 
the contours for kz — where t{kz) has its maximum 
value in (p[). The area between dashed and solid curves 
gives the dispersion in the z-direction of the Fermi sur- 
face in the three-dimensional Brillouin zone. For this 
model band structure, we have evaluated (pl3| ) without 
further approximations (using a three-dimensional nu- 
merical routine over the Brillouin zone as in our previous 
paper [y) and found results (see Fig. 2 in Ref. Q]) for 
the frequency dependence of the c-axis conductivity in 
qualitative agreement with our Fig. |^. 



p,{T) ^ [aduj -^ 0)]-' ^ t-^T 



T' 



4r2 



(40) 



The semiconductor-like upturn behavior at low temper- 
ature is obtained when the second term (interband tran- 
sition) is dominant. 



D. Thermal conductivity 

Based on Wiedemann-Franz law and using ( |37|) a.t ui — 
or directly from (EO), one obtains the c-axis thermal 
conductivity for the plane-plane bilayer 



Kc{T) 



4r2 



(41) 



Once again, Kc{T) will show strong depression at low 
temperature when the interband transition dominates. 





FIG. 6. The structure of the tight-binding bands given 
in Eq. (|42[). The parameters are {ti,i2,/^i, At2,i±} = 
{70, 100, -65, -175, 20} meV and B = 0.45. 
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FIG. 7. Resistivity in the a (top frame, solid line), b (top 
frame, dashed line) and c (bottom frame) directions using 
the tight-binding bands. The a and 6-axis resistivity are for 
t± = 2meV and the c-axis resistivities are for t± — 2, 5, 10, 
and 20 meV from top to bottom. 



In Fig. 1^, we show results for the DC resistivity. The 
top frame refers to the a- and 6-directions, while the lower 
frame is for the c-axis resistivity Pc(T) normalized to its 
value at 300K. Four values oit± have been used for pc{T) 
namely t± — 2, 5, 10 and 20 meV from the top to bottom. 
On comparing this frame with our previous Fig. |3| it is 
clear that there is no qualitative differences between the 
two models and use of the more realistic band structures 
has made only minor difference in our results for the DC 
resistivity. In the top frame of Fig. {t± — 2meV), the 
solid line is for the a-direction and the dashed line for 
the &-direction which is along the chain. While interband 
contributions have been included in the calculations with 
different values of t± , this was found to make only small 
differences on the scale shown and so only a single curve 
is shown. For even higher values of t± , larger differences 
should start showing up. In the case considered, the a- 
b plane resistivity pab{T) is nearly linear in T but not 
quite. We note that our phenomenological assumption 
that the in-plane scattering rate T{T) is linear in T leads 
directly to a linear dependence of Pab{T) vs. T in an 
infinite band model with constant density of states. The 
dashed line in the top frame of Fig. is seen to be close 
to linear but this is not as true for the solid curve (along 
the a-direction) which shows two temperature scale - the 
second presumably related to the Van Hove singularity in 
the two-dimensional tight-binding band structures used 
here. There is also a small interband contribution further 
complicating the situation. 



VII. CONCLUSIONS 

In this paper we have derived approximate expressions 
for the frequency dependence of the c-axis conductivity 
for a bilayer cuprate based on simple free electron bands. 
The simplified results reproduced well our previous nu- 
merical results based on tight-binding bands [||. The 
advantage of the present approach is that analytic re- 
sults can be obtained in several important cases. Besides 
the AC infrared conductivity, we also give results for the 
DC resistivity, the optical conductivity sum rule and the 
thermal conductivity via the Wiedemann-Franz law. In 
the case of the c-axis resistivity, results of full numerical 
calculations based on tight-binding bands not presented 
before are also given and compared with the free elec- 
tron case. In the discussion emphasis is placed on the 
role of the pseudogap which is identified as the minimum 
energy for interband transition possible in a given band 
structure. 
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